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We compute the cross section for e+e~ vvA^ in the general CP-conserving type-II two-Higgs- 
doublet model. We sum the contributions from the "t-channel" e+e~ vuWW — > vuAP graphs and 
"s-channel" e'^e~ ZA" — > uVA'^ graphs, including their interference. Higgs-triangle graphs and 
all box diagrams are included. For many parameter choices, especially those in the decoupling region 
of parameter space (light /i° and mj^o ,mfjo ,mjf± > 2mz) the Higgs-triangle and box diagrams are 
found to be of minor importance, the main contributing loops being the top and bottom quark 
triangle diagrams. The predicted cross section is rather small for tan/3 > 2 and/or m^o > 2mt. 
However, we also show that if parameters are chosen corresponding to large Higgs self-couplings 
then the Higgs-triangle graphs can greatly enhance the cross section. We also demonstrate that 
the SUSY-loop corrections to the bbA° coupling could be such as to greatly enhance this coupling, 
resulting in an enhanced uVA^^ cross section. Complete cross section expressions are given in the 
Appendices. 

PACS numbers: 12.60.Jv, 12.60.Fr, 14.80.Cp, 14.80.Ly 



I. INTRODUCTION 



The Higgs mechanism provides an elegant way to explain electroweak symmetry breaking (EWSB) and the origin 
of the masses of all the observed Standard Model (SM) particles. In many approaches, the symmetry breaking arises 
from a sector involving scalar fields and leaves behind one or more physical Higgs bosons. Detecting and studying 
all such Higgs bosons is one of the major objectives of current and future particle physics experiments. The minimal 
SM contains one SU(2)L-doublet Higgs field, leading to a single physical CP-even Higgs boson after EWSB. However, 
the electroweak scale is not stable with respect to radiative corrections in the minimal SM. Numerous extensions of 
the SM have been proposed to cure this naturalness/hierarchy problem, many of which predict a low-energy effective 
theory with a Higgs sector that contains two (or more) Higgs doublet fields. Our focus here will be on a general 
two-Higgs-doublet model (2HDM) (for a review and references see 0]). In particular, the most promising extension 
of the SM that resolves the naturalness and hierarchy problems is low-energy supersymmctry (SUSY), which must 
contain at least two Higgs doublets. For precise gauge-coupling unification, exactly two doublets arc preferred, as 
incorporated in the minimal supersymmetric Standard Model (MSSM) . The two-doublet Higgs sector of the MSSM 
is predicted to be CP-conserving at tree level and to have type-II fermionic couplings in which one doublet ($i) gives 
mass to down-type quarks and charged leptons while the other doublet ($2) gives mass to up-type quarks. In this 
case there are five physical Higgs bosons: the CP-even /i° and the CP-odd A'^, and the charged pair H^. The 
most important additional parameters of the CP-conserving type-II 2HDM are: (i) tan /? (tan/? = V2/V1 = {^2) 
the ratio of the vacuum expectation values of the neutral members of the two Higgs doublets); and (ii) the mixing 
angle a that diagonalizcs the neutral CP-cven Higgs sector. 

Higgs searches at the CERN LEP II collider have excluded a SM Higgs boson with mass below 114.4 GcV at 95% 
confidence level. In the context of specific choices for the soft-SUSY-breaking parameters at the TeV scale, LEP II can 
be used to exclude a range of MSSM Higgs masses and tan/3. For example, assuming the maximal-top-squark-mixing 
scenario with msusY = 1 TeV, LEP II excludes niho < 91.0 GeV, 771^0 < 91.9 GeV and 0.5 < tan/3 < 2.4 at 95% 
Confidence Level (CL) Searches for top quark decay t H^b at the Fermilab Tevatron exclude tan/3 > 50 

when rrit > mfj± , and searches for the final state bbhP bbbb exclude very high values of tan /3 as a function of the /i° 
mass . Precision electroweak measurements provide only weak constraints on tan /3 . Finally, limits can be placed 
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on tan /? as a function of mfj± based on 7r+, i?+ fx^u decays and on K-K and B-B mixing In the context 
of a type-II model, as defined earlier, these roughly require 1 < tan (3 < 200 for m//± ~ 300 GeV. However, much 
of the MSSM parameter space remains to be explored at future collider experiments. Indeed, for the most general 
MSSM boundary conditions, there is no lower bound on tan/3 from LEP II data. For the most general 2HDM, only 
the presence of two simultaneously light Higgs bosons can be excluded |^ . 

At Run II of the Tevatron, discovery of the light CP-even Higgs boson /i° of the MSSM at the 5(t level is possible 
for m/jo < 120 GeV with 15 fb~^ of integrated luminosity At the CERN Large Hadron Collider (LHC), discovery 
of is virtually guaranteed over all of the MSSM parameter space [Hill El, and measurements of ratios of the the 
partial decay widths in the more prominent decay channels will be possible with precisions on the order of 15% 
A linear e~'"e~ collider could make precision measurements of the couplings with accuracies of a few percent 
In fact, at a linear collider with ^/s > 350 GeV, at least one of the CP-even Higgs bosons of a general 
2HDM (or more complicated Higgs sector) is guaranteed to be detected ^3 in the Zh production mode. In contrast, 
even in the simple 2HDM there is no guarantee that the CP-odd A'^ can be detected, and the situation only worsens 
for more complicated Higgs sectors. Thus, it behooves us to explore every option for A'^ production. 

For a CP-conserving Higgs sector, production of a single A'^ via loop-induced processes could prove critical to a full 
exploration of the Higgs sector. This is because the most useful tree-level mechanisms for single Higgs production rely 
on a substantial Higgs coupling to ZZ or WW pairs. Such couplings are absent at tree-level for the purely CP-odd A'^. 
If the Higgs sector is CP-violating, then the neutral Higgs bosons will mix with one another and, in general, all will have 
substantial tree-level ZZ and WW couplings. As a result, in the case of a CP-violating Higgs sector, loop-induced A'^ 
production mechanisms would probably not be very important. Thus, loop- induced A^ production is mainly of interest 
for a CP-conserving Higgs sector. We note that there are substantial phenomenological reasons for believing that the 
Higgs sector will prove to be CP-conserving. In particular, CP-conservation is the most straightforward approach to 
avoiding conflict with the constraints coming from the anomalous magnetic moment of the muon 

(5m ~ 2) [ig( and the 

non-observation of electric dipole moments (EDMs). Still, we must note that even in the MSSM context substantial 
CP violation could be introduced at the loop level if the soft-SUSY-breaking parameters have phases, and that a 
CP-violating Higgs sector can be consistent with the EDM and {g^ — 2) constraints if there are carefully orchestrated 
cancellations between CP-violating contributions to these observables |23| . 

Let us review in detail the difficulties associated with producing and detecting a purely CP-odd A'^. At a hadron 
collider, the absence of tree-level ZZA'^ and WWA'^ couplings implies that: (i) the W* W production mode is 
suppressed, as particularly relevant at the Tevatron; (ii) the — > WW and gg Al^ ZZ production/decay 
modes (the "gold-plated" processes for discovery of a heavy SM-like Higgs boson) have very low rates because the 
branching fractions B{AP ZZ) and B{A'^ WW) are small; and (iii) at the LHC, the gg ^ A'^ ^ 77 rate is 
numerically small. With regard to the latter, we note that had r(^° 77) been of SM-like size, the absence of 
tree-level A'^ WW,ZZ decays would have implied substantial B{A^ — > 77) and a useful gg A'^ 77 rate 
even for large m^o. However, the absence of the M^-loop (the largest contribution in the case of a SM-like Higgs) 
results in an even greater suppression of T{A'^ 77) than of F^^. At an e+e^ collider, the e+e^ Z* ^ ZAl^ and 
e+e" vijW*W* — *■ vvJ^ processes are only present at the one-loop level. 

In general, the A^' can be pair produced at tree-level. However, the rates for Higgs pair production are generally 
too small for observation at the Tevatron and LHC since they are electroweak in strength and must compete against 
enormous QCD backgrounds. Pair production is, however, potentially useful at a e^e^ machine. Such processes 
include e+e" ^ Z* ^ ZA°A° IHlIl, e+e" ^ wW*W* uVA^A^ [Hill and e+e" -> Z* -> H^A^ or h'^A^ 
(see Q). However, all of these processes can be simultaneously suppressed by kinematics and/or small couphngs. In 
particular, this occurs in the decoupling limit of a 2HDM that typically arises when m^o > 2mz 131. In the decoupling 
limit m^o ^ rn^o ^ miio, implying that ^/s > 2m^o is required for H'^A^, ZA'^A^ and i/VA'^A^ production (all of 
which would otherwise have large cross sections since the ZZA^A^ and WWA^A^ couplings arc fixed by the standard 
quadratic gauge couplings of the A'^ and the ZH^A^ coupling is maximal in this limit), while h^A^ production is 
strongly suppressed in the decoupling limit by a factor of order ra%/m\o in the square of the hPA^Z coupling. This 
decoupling limit is automatic in the context of the MSSM and is quite natural in the case of a more general 2HDM. 
In the general 2HDM there are other scenarios not related to this standard decoupling limit in which detection of the 
on its own would also be critical. In particular, it is possible to choose Higgs sector parameters in such a way 
that the A^ is the only light Higgs boson while maintaining consistency with precision electroweak measurements |25j 
(see also |23|). Further, it is possible that a light A^ could explain part of the observed discrepancy between the SM 
prediction for (g^ — 2) and the experimentally measured value p^ . For all the above reasons, it is important to assess 
more carefully the various possible mechanisms for single A^ production. 

Consider first the tree-level processes for single A^ production. At both the LHC and at an e+e~ collider, the 
only relevant tree-level processes for single A^ production are tiA^ and bbA^ production. At the LHC, it has long 
been established that these do not yield an observable signal in a wedge-shaped region of the (m^o , tan (3) parameter 
space. This wedge spans a range of moderate tan (3 values for m^o > 200 GeV and becomes increasingly broad as 
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forbidden) there is no known means for detecting the A'^ using tree-level production mechanisms when tan /? is not 
large enough for the bhA'^ process to have an observable rate. In this situation, we must turn to loop-induced A° 
production mechanisms. 

One possibility is to build a photon collider at the e+e~ collider and look for 77 — > A'^ via t, b and charged Higgs 
loops Ji, .32_, In particular, the recent realistic study of Ref. [s^ found that three years of running at a 630 GeV 
e+e" LC in the photon collider mode could provide a 4(t signal for a significant fraction of the (m^o , tan (3) LHC wedge 
region with m^o < 500 GeV. If m^o could be rou ghly gue ssed, e.g., based on the precise measurements of deviations 
of the couplings from their SM values [T^ l3^ I35ll3a . l37l | obtainable at a linear e+e" collider, then the energy of the 
photon collider could be chosen to optimize the A" production cross section resulting in a faster discovery. One-loop 
production possibilities in e+e" collisions include e+e~ — *■ jA^, e+e" ^ Z* — > ZA'^ and e+e~ — *■ I'VA'^ via WW 
fusion. (The latter two interfere since Z vv.) The e+e^ 7^4° process has been computed in Refs. j^MES- Results 
for the ZA^ final state (computed assuming the absence of supersymmetric particle loops) appear in Refs. [2^ Is^ . 
The possible enhancement of the ZA^ rate when SUSY particles are present is discussed in [i^l . 

In this paper, we give a complete calculation of the W boson fusion process e+e^ vVW*W* vvA^ in the 
general CP-conserving 2HDM, which first occurs at one-loop since there is no tree level W~^W~ A^ coupling. We 
include the process e^e" ^ Z* —t ZA^ vV^ (first computed in Ref. H^), which leads to the same final state 
and thus interferes with the W boson fusion process. The process e+e~ vvJ^ has also been computed recently in 
Ref. After a review of the structure of the general 2HDM in Sec. in] we present the relevant Feynman diagrams 
and formulae involved in our calculations in Sec. IIIII We present numerical results in Sec. IIVI 

In Sec. lIV AI we employ the tree- level MSSM two-doublet sector (see 0) as a benchmark for our study. The Higgs 
sector of the MSSM is a type-II 2IIDM in which the quartic couplings of the two Higgs doublet fields are fixed at tree 
level by the gauge couplings. In this case, all Higgs sector parameters are determined by just the two parameters m^o 
and tan/3, and for mj^a > mz the 2HDM quickly approaches the decoupling limit, described earlier, in which all Higgs 
self-couplings are small. We compare the full 2HDM results including the loop contributions from quarks, Higgs and 
gauge bosons with those obtained by including only the top and bottom quark loops, as computed in |22l |. We show 
that the e+e" vvJ^ process could provide a viable signal for tan/3 < 1, thus covering part of the region where 
discovery using tree-level processes is not possible. Further, for such tan /3 values the vv^ rate would provide 
a very sensitive measurement of tan /3. In contrast to the above, the vVJ^ cross section is typically quite small for 
tan /? > 2 for parameter choices based on the tree-level MSSM Higgs sector potential. Of course, in the full MSSM, 
radiative corrections to the tree-level masses should be incorporated as should the contributions with superparticles 
running in the loop. However, the sparticle loops are in general suppressed by the heavy sparticle masses. A full study 
of the MSSM, including all the superparticle loop contributions and radiative corrections to the Higgs potential, will 
appear in 42]. 

Sec. IIV Bl focuses on the general 2HDM with Higgs potential parameters outside the decoupling regime. We find 
that if the Higgs self-couplings are large (as possible when 2HDM parameters are chosen to lie in a non-decoupling 
regime) then the vvAl' cross section can be greatly enhanced by Higgs triangle diagrams even for such large /i", 
masses that none of these latter Higgs bosons could be directly observed. Indeed, for lower -01^0 values (m^o < 2mt), 
the rate is sufficiently enhanced when tan/3 > 10 that the vvAS' events could have a detectable rate and provide direct 
evidence for the large Higgs self-couplings. This probe of the Higgs self-couplings would be especially powerful if some 
of the other Higgs bosons have themselves been directly observed. 

In Sec. IIV CI we illustrate one unusual possibility in the full MSSM context, namely that the hbJ^ coupling could 
be greatly enhanced by non-decoupling SUSY particle loops, resulting in a huge enhancement for the e^e~ — > vvAP 



Finally, Sec. is reserved for our conclusions. In the Appendices, we collect the complete matrix elements for the 
various Feynman diagram contributions. 



V = m\^^\^i + m\^^\<^2 - K2$t$2 + h.c] + iAi($t$i)2 + \\2{^\^2f + A3($I$l)($^$2) + A4(*I$2)($J$l) 



cross section. 



II. THE CP-CONSERVING 2HDM 



We adopt the conventions of 01 for the 2HDM. Let $1 and $2 denote two complex SU(2)i-doublet scalar fields 
with hypercharge Y ^ \. The most general gauge invariant scalar potential is given by 




(1) 
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The terms proportional to Xq and A7 lead to flavor-changing neutral current interactions (FCNCs) and will be set 
to zero. This can be achieved by imposing a discrete symmetry <i>i — > — $1 on V. However, we allow for a soft 
(dimension- two) breaking of this symmetry through 771^2 7^ 0.^ If Ag = A7 = but mfj 7^ 0, the soft breaking of the 
discrete symmetry generates finite Higgs-mediated FCNCs at one loop. The tree-level supersymmetric form of V is 
obtained from Eq. (QJ by the substitutions: 

Ai=A2 = i(.g2+.g'2), A3 = i(,g2-<7"), ^4^-^^ X5 = Xe = X7 ^ , (2) 

where g and g' arc the SU(2)^ and U(l)y gauge couplings, respectively. In general, ^^'^ '^5 (and, if present, Ag 
and A7) can be complex. However, we explicitly exclude such CP-violating phases by choosing all coefficients in Eq. 
to be real and such that spontaneous CP violation is absent. For details, see Ref. |23|. The scalar fields will develop 
non-zero vacuum expectation values if the mass matrix has at least one negative eigenvalue. Imposing CP invariance 
and U(1)em gauge symmetry, the minimum of the potential corresponds to the following vacuum expectation values: 

where Vi are real in the absence of explicit and spontaneous CP violation. The minimization conditions on the 
potential can then be used to determine to^j and 77132 in terms of the other parameters (with Ag = Ay = 0): 



mj^ = ^12^/3 - [Xicjj + X345S}] , m22 = ml2ti3^ ~ [A2s| + A345CI] , (4) 
where we have defined: 

A345 = A3 + A4 + A5 , i/3 = tan /3 = — , (5) 

and = v\ + V2 ~ imf^/g'^ ~ (246 GcV)'^ . It is always possible to choose the phases of the Higgs doublet fields 
such that both vi and V2 arc positive, implying that we can take < (3 < Tr/2. With Ag = A7 = 0, all but one of the 
eight free parameters in Eq. can be fixed after EWSB in terms of v, tan/3, the four physical Higgs masses, and 
the mixing angle a required to diagonalize the neutral CP-even Higgs sector. In our numerical analysis below, we use 
the parameter set m^io, ttiho, to^o, mjj±, a, tan/3 and A5. (v is of course fixed by the measured values of mw and 
g.) The relations (D.20)-(D.23) of ^ with Ag = A7 = then give the other X^ as: 

_ mjjoc'^ + mlos'^ - m\os'p „ 
Ai - -5-2 Agt^; [b) 

2 2 2 2 2 

A2 - -T-2 ^5t^ ; (7) 

_ {m'\ja-m\o)SaCa + {2m^^±-m\o)s(iCi3 



(8) 



A. ^ 2(7.^0 -m^,) ^^^^ 

The mass parameters of the Higgs potential are given by (D.17), (D.24) and (D.25) of [2^ : 

™12 = spCisiX^v"^ +m\o)] (10) 
mil = ~-^{rn]joCaCp-a-m\oSaSp-a) + s'liX^v'^ + m\o)] (11) 

77122 = --^^{m\oCaSp-a+rn]joSaCl^-c) + c^piX^v"^ +m\o) ■ (12) 

In the supersymmetric limit, the A^ are determined in terms of gauge couplings as given in Eq. ((SJ and the Higgs 
sector is then entirely specified at tree-level by the two free parameters tan (3 and 777,^0 . 



^ This discrete symmetry is also employed to restrict the Higgs-fcrmion couplings so that no tree-level Higgs-mediated FCNCs are present. 
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FIG. 1: Feynman diagrams for 2HDM contributions to e vvAP via f-channel processes. Here Si denotes Higgs and 
Goldstone bosons, Vi denotes gauge bosons. 

There are two types of 2HDM, depending on which Higgs field is responsible for the masses of quarks and leptons. 
In the type-I 2HDM, $2 gives masses to both quarks and leptons. In the type-II 2IIDM, $2 couples to up-type quarks 
while $1 couples to both down-type quarks and charged leptons. Consequently, the Yukawa couplings of the quarks 
and leptons to the Higgs bosons are different in these two cases. For the ^ we find L =- ii/ffj^fA'^ where: 



yt = — cot p Ub = — cot p ijr ~ cot p Type I. 

V V ' V 



yt = — - cot f3 Vb = — - tan /3 I/t — — - tan /3 Type II. 

V V V 



nib 



(13) 
(14) 



The MSSM is required to have type-II couplings and our analysis will also assume type-II couplings for the general 
2HDM case. 



III. FORMALISM 



In our analysis, we adopt the renormalization scheme of Ref. The diagrams contributing to e+e" — *■ vVA^ in 
the 2HDM via f-channel processes are shown in Fig.^ where we have neglected all the diagrams that are proportional 
to the small electron Yukawa couplings. Diagrams (l)-(7) are the (finite) triangle loop corrections to the effective 
VF+Vl^"^^ coupling, and diagrams (8)-(ll) are the box diagrams. For on-shell W bosons, the sum of diagrams (3)-(7) 
with Higgs bosons and vector bosons running in the loop must exactly cancel li^j. But, in the present context, the W 
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+ crossed 
diagram 



Z 



(3) (4) 

FIG. 2: Feynman diagrams for 2HDM contributions to e~ vvJ^ via s-channel processes. 

bosons are virtual and the sum of these diagrams is non-zero. There is no WW counterterm contribution since the 
WW vertex is finite. In addition, there are no tadpole contributions since we have set the renormalized tadpoles 
to be zero. We shall refer to this collection of diagrams as the t-channel diagrams. 

There are additional diagrams related to s-channel 7, Z exchange; these are shown in Fig. |21 Diagrams (1) and 
(2) are the triangle loop contributions to the and ZZA couplings (which are zero at tree-level), followed by 

Z —f vV. Diagrams (3) and (4) are related box diagrams; photon exchanges do not appear in these diagrams because 
of the absence of ^vv and ^ZS vertices. If the final exchanged Z connected to vv is on-shell in diagrams (l)-(3), then 
the calculation is equivalent to e+e^ ZJ^. This process has been calculated for the 2IIDM and the full MSSM 
in [39I l40l | and the production cross section found to be small. The Z resonant contribution to the vVJ^ final state 
can be separated experimentally by detecting the A" in a visible final state (e.g. 66), reconstructing the mass of the 
vV recoiling opposite the A° and removing events in which the reconstructed mass is near mz- However, far off-shell 
intermediate Z bosons can potentially give s-channel contributions that interfere with the ^-channel contributions to 
the vvAl' final state, and must be included in our calculation. 

One possible concern is that including the Z decay width in the Z propagators in the diagrams of Fig. |21 might 
spoil the gauge invariance of the calculation. We checked explicitly and found that this is not the case. In particular, 
the sum of diagrams (1) and (2) in Fig. |21is gauge invariant on its own. Also, there is a cancellation between the 
gauge-dependent part in diagram (3) and the corresponding crossed diagram. A similar cancellation occurs between 
the gauge-dependent part of diagram (4) and its crossed counterpart. All these cancellations occur as a result of 
numerator algebra and are independent of the Z width appearing in the Z propagators. 

Another check of the gauge independence of our fixed width scheme for the Z boson is to compare the numerical 
results to those of the "factorization scheme", which is guaranteed to be gauge independent. Following, e.g., the 
discussion in Ref. |45l |. the one-loop matrix element in the factorization scheme is given by^ 

M = Mr.^o. (15) 



^ For processes that are nonzero at tree level, care must be taken to avoid double-counting the Z width. This is not a concern here since 
the tree level matrix element is zero. 
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(a) (6) (c) (d) (e) 




FIG. 3: Feynman diagrams for the contribution to e'^e — > uuA'^ through Z-A° and j-A^ mixing self-energies. The blob 
denotes the renormalized Z-A^' and 'y-A^^ mixing. 

The factorization scheme sets the nonresonant diagram (4) of Fig. [3 to zero when s = A/|-, which leads to an effect 
of 0{arz /Mz) ^ 0{a^); since this is a higher-order effect, it should be small. Om' results in the fixed width scheme 
agree numerically with those of the factorization scheme to within the precision of our phase space integration. 

Another set of diagrams involving Z — Al' mixing is given in Fig. O The contribution from diagram (a), in which 
the W~^W~ couple to a virtual Z* which then mixes with the A" via (infinite) one-loop diagrams, is zero for on-shell 
W bosons With virtual W bosons as in our case, diagram (a) is non-zero. However, it is not gauge invariant on 
its own. Additional diagrams (b)-(e) have to be taken into account. The sum of all these diagrams gives zero as a 
consequence of gauge invariance. This can be seen as follows. The one-particlc-irrcducible (IPI) two-point function 
for Z — mixing is defined as —ip^Y^z-AoiPA), where is also the off-shell Z momentum. Gauge invariance for 
the Z tells us that after summing over all possible Z attachments we must have p^o^p = 0, where Ai_i is the full 
amplitude that would be dotted into the Z propagator. In some renormalization schemes, one can also have 7 — ^4*^ 
mixing, in which case a similar argument guarantees that after a complete sum over the diagrams (a), (c) and (e) 
with a virtual 'y*OTie gets zero, just as in the Z* case. Of course, there are renormalization schemes (such as that 
chosen in Ref. |43|) in which there is no 7 — mixing and this issue does not arise. Similarly, the sum of s-channel 
diagrams (f)-(i) also gives zero. 

The diagrams shown in Fig. ^ do not contribute to our calculation. For diagrams (a) and (b), the one-loop 
— > H'^ mixing graph must be proportional to = ~ P^- (for example) , which gives zero when acting on 
the V^iz77^e vertex by the equation of motion. (Here we take the approximation that both electron and neutrino are 
massless.) Similarly, diagrams (c) and (d) also vanish. 

Let us define p^- , Pe+, p^, pjj, ki = Pu — Pe- k2 — Ptj — Pe+ to be the momentum for the incoming electron, 
positron, outgoing neutrino, anti-neutrino, intermediate W'^ and W~ ^ respectively. The i-channel matrix clement for 
e+e" vvA^ can be decomposed as follows (using ^ = 1 Feynman- 't Hooft gauge): 

_ v{Pe+ , Se+hi^PLvjPU, S-) u{p^, S^)jf,PLu{p^- , S^- ) 
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(a) (b) (c) id) 

FIG. 4: Contributions to e'^'e" vVJ^ tiiat we neglect. The blob denotes the renormalized W^-U^ mixing in diagrams (a) 
and (b) and the renormalized {j, Z)-{h° , mixing in diagrams (c) and (d). 

X [Fe'^^^'-h afc2 + Gg^^^ + H.p'^^p^^ + H^p'^p^ + H.^-p'^^p^) + Hii-pX^ )] ' (16) 

where = (1 — 75)/2. For the triangle loops [Fig. diagrams (l)-(7)], only the combination kik!^ appears for the 
Hi term. Therefore, the terms inside the square brackets in Eq. H16() can be simplified as 

[Fe^'^'-'h + + Hk'(kip^ , (17) 

where H = Hi = H2 ~ = H4. Recall that the index /i (i/) is associated with the going into the 

(e"*") line. Analytical formulae for each Fcynman diagram contribution to F, G, Hi are summarized in Appendix IbI 
In the on-shcU limit where kf = fc| = 'm'w^ G and H arc both zero 01; which can be seen explicitly in the analytical 
formulas. 

The s-channel matrix element can be written as two non-intcrfcring pieces: Ais = -Msie^s^) + s{e^e~j^) , where 

i even i odd 

The definitions of the operators Oi and the contribution of Fig.|21to A4i are given in Appendix IHl 
The spin averaged matrix element squared is: 

lE . I-^P - lE . 1-^*1' + !E . \Msiele^,)\' + ji: . \Msie],et)\' 

4 spiJis 4 — ' spins 4 ^"^^ spins 4 ^"^^ spins 

+ lE ■ MtMs{ele+y (19) 

4 ^ — ^ spins ^ it 

which is used in the cross section calculations. The factor of 3 in the second and third terms represents the sum over 
the three neutrino flavors in the s-channel contribution. Notice that only the matrix element A^s(e^eJ) interferes 
with the i-channel diagrams, while the other s-channcl matrix element does not. The explicit expression for the pieces 
in Eq. (|19|l is given in Appendix IdI 

Our numerical computations were performed using the LoopTools package j4fij . We thus write the Appendices 
using the notation of LoopTools for the one-loop integrals. 

IV. NUMERICAL RESULTS 

A. 2HDM with tree-level MSSM mass and coupling relations 

In this section, wc give results for the e'^e^ ^ vVA^ cross section as a function of m^o and tan (3 assuming a typc-II 
2IIDM with the tree-level MSSM constraints on the quartic couplings in the Higgs potential. As a result, the masses 
771^0 , rrifjo and 77i^± and the mixing angle a of the CP-even sector are all fixed in terms of m^o and tan (3 by the 
tree-level relations of the MSSM. The tree-level MSSM couplings lead to a theoretical upper bound on the mass of the 
lighter CP-even of m/jo < mz (which is increased to ^ 135 GeV by radiative corrections 01)- decoupling 
region of large m^o (typically m^o > 2mz is large enough), the only light Higgs boson is the CP-even whose 
couplings to the SM particles approach their SM values. The other Higgs bosons and H^ can be as heavy as 

a TeV. The heavy Higgs bosons are nearly degenerate in mass, with mass splittings of the order of m^/m^o. 
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Ecm=2mA (GeV) 

FIG. 5: Behavior of the t- and s-channel contributions to e'^e~ vVA:^ with Ecm- Here tan/? — 2.5 and m^o — i5cm/2. The 
line labeled "total" shows the sum of the t- and s-channel contributions including their interference. The s-channel results 
including gauge/Higgs contributions are indistinguishable from the s-channel results with only top/bottom loops. 

Our choice of the tree-level MSSM Higgs sector can be thought of simply as a representative model choice within 
the 2HDM that gives decoupling and a SM-like ft," as m^o gets large. It will allow us to explain general features 
of the cross section and how they depend upon mj^a and tan /3. We will only examine results for tan /? > 1 in this 
section. Our focus will be on situations in which m^o > -Ecm/S, implying that pair production of the (e.g., 
e+e" — > vv^J^^ will be kinematically forbidden, as will e+e" A^iJ", since ~ m^o for MSSM-like mass 

relations. The e^e~ hb^ process will also be strongly suppressed for m^o > Iniz- Thus, we are considering 
situations in which ^ discovery might only be possible via the (one-loop) single production mechanisms that we 
consider. 

The first feature of interest is that the s-channel and i-channel contributions to the cross section have different 
behavior as the collider center-of-mass energy increases, as shown in Fig. [31 In particular, for m^o = £'cm/2, the 
s-channel contribution dominates for E^m ^ 650 GeV, while the t-cliannel contribution dominates for Ecm Si 700 
GeV. To some extent the s- and i-channel contributions can be separated experimentally. For example, the s-channel 
contributions can be isolated by looking at ZA'^ final states in which the Z decays to an observable final state, such 
as qq or £'^£~ . The i-channel contributions can be isolated to a large extent by looking at the A'' in a visible final 
state decay mode (for example, bb or r"'"r~), reconstructing the mass recoiling against the A'^, and demanding that 
this recoil mass not be close to mz ■ In this latter case, this selection procedure would reduce somewhat the <-channcl 
cross sections presented (which are integrated over all recoil masses). 

Perhaps most importantly, we find that the cross sections are quite small, generally below 0.001 fb. This holds 
even for rather low values of m^o . The cross section as a function of m^o is shown in Fig. This figure shows the 
expected peak in a(vvA^) at to^o ~ 2mt from the top quarks in the loop going on shell, followed by a rapid fall as one 
approaches the kinematic limit at m/^a = -^s. Notice that for Ecm — 500 GeV, the s-channcl contribution dominates, 
while for Ecm = 1000 GeV, the t-channel contribution dominates, as expected from Fig. [S] Similar results were 
presented in Ref. |4ll | . While we include all three flavors of neutrinos in the final state since they are experimentally 
indistinguishable, Ref. included only VcVe in the final state, leading to an s-channel cross section smaller by a 
factor of 3 than our result. Taking this into account, our results are in rough agreement with those of Ref. For 
example, for the point Ecm = 500 GeV, vtla = 100 GeV, tan [3 = 2.5, and one neutrino flavor, our result is about a 
factor of two larger than that shown in Fig. 4 of Ref. Part of the discrepancy, a factor of (137/128)"* = 1.3, is 
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FIG. 6: Cross section for e+e^ — > uVA^^ as a function of m^o, for tan/3 — 2.5 and Ecm. ~ 500 GeV (left) and 1000 GeV (right). 
The hne labeled "total" shows the sum of the t- and s-channel contributions including their interference. As in the previous 
figure, the s-channel top/bottom only and top/bottom+gauge/Higgs lines are essentially identical. In the Ecm = 1000 GeV 
plot, we note that the "f-channel" lines (whether top/bottom only or top/bottom+gauge/Higgs) are always below the "total" 
lines in the m^o ~ 200 GeV region. Thus, for instance, the upper dashed line is the total "top/bottom only" result while the 
lower dashed line is the t-channel "top/bottom only" result. 



explained by the use of a = 1/137 in Ref. versus our use of a ~ 1/128.'^ 

The dependence of the cross section on tan/3 is shown in Fig. [7| This plot clearly shows that for the MSSM-like 




0.2 0.5 1.0 2.0 5.0 10.0 20.0 50.0 0.2 0.5 1.0 2.0 5.0 10.0 20.0 50.0 

tan/S tan/S 



FIG. 7: tan/3 dependence of the cross section for e'^e~ ui/A°, for Ecm = 2m^o = 500 GeV (left) and 1000 GeV (right). 
As in the previous two figures, the s-channel curves with and without the gauge/Higgs contributions are indistinguishable. 
In the Ecm = 500 GeV {Ecm = 1000 GeV) case, the "total" curves are slightly higher than the "s-ch." ("t-ch.") curves 
for both "top/bottom only" and for "top/bottom+gauge/Higgs". The "top/bottom only" curves are slightly above the 
"top/bottom+gauge/Higgs" curves in both the "s-ch." ("t-ch.") and "total" cases. In this figure, we have used the shorthand 
"ch." for "channel". 



parameter relations, detection of the lyVA'^ final state will only be possible if tan/3 < 1. The cross section falls like 
a power law with increasing tan /3. This is due to the fcrmion triangle diagrams with a top quark coupling to 
(diagram (1) in Fig.Hand the t-loop case for diagrams (1) and (2) in Fig.|21), which dominate at low tan/3 and give 



^ This was pointed out in a private communication from the author of |41| in which he also gives numbers (that include the a correction) 
that are in close agreement with ours. 



11 



a cross section proportional to ~ (tan/3)~^. At large values of tan/3 > 30, the fermion triangle diagrams with a 
bottom quark couphng to (diagram (2) in Fig.^and the 6-loop case for diagrams (1) and (2) in Fig.|2Jl begin to 
contribute significantly and affect the dependence on tan /3, since these diagrams give a cross section proportional to 
2/2^(tan/?)2. 



e*e ^7A: Vi"=500 GeV e*e -*vv k: Vs"=500 GeV e*e -^yA: Vi'=800 GeV e*e -*vv k: Vi'=800 GeV 
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FIG. 8; The e^e~ — > 7^4" and vVA^^ cross sections as a function of m^o for y/s = _Ecm ~ 500 GeV and 800 GeV, for 
tan/3 — 0.5,1,5,20,50. We employ the tree-level MSSM parameterization. The e^e" vVA^ cross section includes all 
contributions. 



Another perspective on these results, and a comparison to the 
For the tree-level MSSM type of 2HDM parameter choices, the e+e" 



^AP process HE US is presented in Fig. |H1 
^A^ process would probably lead to earlier 



discovery of the A^ than would the e^e —> vvA^ assuming both have small background. However, as described 
in the next section, the e+e^ vVA^ can be enhanced for non-decoupling 2HDM parameter choices that lead to 



large Higgs self-couplings, whereas the e+e^ 7^° cross section is not sensitive to Higgs self-couplings and would 
not be enhanced in such a parameter regime. Note also that the advantage of the e^e~ ^A^ process over the 
e+e^ — > vvA^ process decreases slowly with increasing -Ecm, as seen from the figure by comparing the results for 



-Br. 



500 GeV to those for & 



800 GeV. 
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FIG. 9: Contribution of the gauge and Higgs boson loops relative to the top and bottom quark loops, as 
Plotted is the gauge and Higgs contribution normalized to the top and bottom quark contribution, {atotai - 



I function of m^o . 
iytb)/o-tb- 



The relative size of the gauge and Higgs boson contributions compared to the top and bottom quark contributions 
is illustrated in Fig. |51 For tan/3 = 1, the gauge and Higgs boson loops vanish. For larger values of tan/3, the 
contribution of the gauge and Higgs boson loops relative to that of the top and bottom quark loops typically increases 
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with increasing tan/3. At Ecm = 500 GeV, the gauge and Higgs boson contributions can be quite significant at low 
m^o < 100 GeV, especially for larger values of tan f3. However, in the MSSM, m^o < 92 GeV is excluded by the LEP 
II data 0, Q ■ For larger values of m^o , the gauge and Higgs loops interfere destructively with the dominant top and 
bottom quark loops, resulting in a reduction of the cross section by less than 10%. At Ecm = 1000 GeV, in contrast, 
the destructive interference of the gauge and Higgs boson loops with the top and bottom quark loops is much more 
significant for mj^o below the top quark pair production threshold of 350 GeV, suppressing the cross section by as 
much as 30% for m^o ~ 200 GeV. For m^o > 350 GeV, the gauge and Higgs boson loops reduce the cross section 
by less than 10%. The change in the relative size of the gauge and Higgs boson loops at different center-of-mass 
energies can be understood as follows. At Ecm = 500 GeV, the cross section is dominated by the s-channel diagrams. 
The gauge and Higgs boson contributions to the s-channel matrix element come only from box diagrams [diagrams 
(3) and (4) in Fig. At Ecm = 1000 GeV, the cross section is dominated by the t-channel diagrams. The gauge 
and Higgs boson contributions to the i-channel matrix element come from both triangle diagrams and box diagrams 
[diagrams (3)-(7) and (8)-(ll), respectively, in Fig. The gauge and Higgs boson triangle diagrams in general give 
larger contributions to the cross section than the box diagrams, leading to larger gauge and Higgs boson contributions 
to the i-channel process than to the s-channel process. This behavior can also be seen in Fig. IHl 

B. General 2HDM 

There is, of course, much more freedom in the general 2HDM than we have allowed for in the previous section 
First, there is the possibility of allowing for type-I fermionic couplings as opposed to the typc-II fermionic couplings 
employed so far. As shown in Eq. H13|) . the type-I ti coupling is the same as the type-II coupling. This implies that 
the dominant t loop contribution to the cross section will be unchanged. The type-I bb coupling is proportional 

to cot f3 as opposed to tan j3 for type-II coupling; this means that the &-loop contribution to the i/VA'^ cross section 
is never important for type-I couplings. Numerically, this implies that the leveling off of a{in/A'^) at tan/3 > 40 in 
Fig. (and eventual rise at still larger tan/3) would not take place for type-I couplings. There is also the possibility 
of so-called type-Ill fermion couplings in which the two Higgs doublets both couple to up and down type quarks. In 
general, such couplings yield flavor changing neutral currents that are too large compared to existing experimental 
constraints. In addition, the numerical modifications to the type-II predictions already given would not be large. 
Thus, we do not consider typc-I or typc-III couplings further. 

A second variation in the general 2HDM context is to allow for CP-violating couplings. As explained in the 
introduction, we have chosen not to explore this possibility here as it makes the A° less unique and because considerable 
cancellations between CP- violating contributions deriving from the Higgs sector are required for consistency between 
the computed EDMs and (g^ — 2) and experimental data. If the Higgs sector is CP-violating, all the neutral Higgs 
bosons mix and will all have some level of VV coupling. In most scenarios, all three of the neutral Higgs bosons hi 
[i = 1,2,3) would be easily detected in Zhi production and the tree- level contributions to the vvhi {i = 1,2,3) cross 
sections would all be considerably larger than the one-loop contributions . 

The most interesting issue in the general CP-conserving 2HDM context is the extent to which the Higgs self- 
couplings could deviate from those in our previous MSSM-like analysis, so that (j{i'VA'^) might be substantially 
increased or decreased by the Higgs boson loops compared to the value obtained from the t and b loops. In Fig.|51we 
found that with MSSM-like couplings, the gauge and Higgs boson loops could change the cross section by as much as 
70% compared to the t and b quark contributions for low values of m^o, or by up to 30% for m^o > mz- Here, we 
explore the effect of removing the MSSM constraint on the Higgs boson self-couplings Ai, while still requiring that 
they remain perturbative. (Following [23|, we define perturbativity by the requirement Ai/(47r) < 0(1).) Another 
constraint on the 2HDM parameters derives from precision electroweak data, as conveniently summarized by the S 
and T parameters. We will explore the extent to which the vvA^ cross section can be enhanced via the triangle 
graphs involving Higgs self-couplings while remaining consistent with the perturbativity and S, T constraints. In the 
preceding section, we considered parameter choices that correspond to rapid decoupling as m^o increases beyond 
0{mz)- The Higgs self-coupling effects are likely to be most significant in those regions of parameter space that are 
far from the decoupling limit. 

The triangle diagrams involving one or more internal gauge bosons, Fig.^diagrams (4)-(7), are controlled by VSS 
and VV S couplings which are determined by gauge invariance. The size of these couplings is limited, as shown in 
Tables HTl and HTll Any significant enhancement must come from the purely Higgs loop graph of Fig. ^ diagram (3), 
which is controlled by the Higgs self-couplings given in Table n For the dia grams in Fig. ^ diagram (3) involving one 
or more Goldstonc bosons, the couplings arc determined by gauge invariance purely in terms of the Higgs masses, 
the angle combination (3 — a and the weak mixing angle Ow (see Table IJ). For the diagrams with h^H^A^ or 
H'^H^A^ running in the loop, on the other hand, the couplings gAOA°h° ^tnd gAOA^H" ^-re determined by the invariant 
combinations of the Aj, denoted by At and \u in Ref. 0, given in Eqs. (jB13|) and (jB14|l . These couplings are 
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free to vary in the general 2HDM once the MSSM constraints are removed. We will employ the ratios Xjj / (An) and 
Xt/{Att) to quantify the perturbativity of the self-couplings. Numbers much larger than 1 for these ratios imply that 
the self-couplings are becoming non-perturbative and that higher order corrections to the results we obtain could be 
large. 

Our procedure will be to choose values for the Higgs masses m^jo , m^o and mj^o and for tan /3, and then scan over 
a and m/^i selecting points that are consistent with the experimental S, T values at the 95% CL. As in Ref. (see 
also HSl), for simplicity we will restrict A4 relative to A5 by requiring A4 = — A5, implying = (m^± — m\o)- 

This choice makes it relatively easy to find values for the other parameters that give good agreement with precision 
electroweak 5, T data. Indeed, as discussed in Ref. p5| . even if we choose large values for m/jo and m/^o (in particular, 
beyond the kinematic reach of the linear collider), it is nonetheless possible to choose a and mfj± in such a way that 
the iS*, T values are within the 95% CL ellipse. The key is to have > m^o, m^/o by a small amount (typically 10 

to 30 GeV) in such a way that the large negative AT generated by the heavy neutral scalar(s) with substantial VV 
coupling is compensated by an even larger positive AT contribution coming from the 'mjj± — m^o and/or m//± — m^o 
mass difference. As one varies a at fixed tan/3, it is generally possible to adjust the value of 'mjj± in such a way as 
to remain within the (upper right hand segment of the) 95% CL S-T ellipse. 

We show the results of this procedure in Fig. ^| We compute aiyvJ^) for m^o = 100, 350, and 600 GeV as a 
function of tan/3, assuming a collider energy of -y/s = 800 GeV. We take to/jO = 800 GeV and m^o ~ 900 GeV and 
scan over < a < tt in 10 steps and over to//± in steps of 10 GeV beginning with m/fo as the lowest value. In 
particular, all the Higgs masses other than m^o are chosen such that none of the other Higgs bosons can be produced 
for the assumed collider energy. We plot only those points for which the S*, T values are within the 95% CL precision 
electroweak ellipse. In the left-hand plots of Fig. ^| we give a(vvJ^) with and without including the Feynman 
diagrams containing Higgs self-couplings. In the right-hand plots of Fig. ^| we plot the corresponding values of 
AT/(47r) and A£//(47r) (without attempting a point-by-point identification). We observe that the Higgs self-coupling 
diagrams can have a very large effect on the cross section at large tan /3 if one is willing to accept values of \t,u / (47r) of 
order 2 to 3. At large tan/3, the cross section is very substantially enhanced by the self-coupling graphs and might be 
visible with L = 1000 fb~^ of integrated luminosity. Certainly, detection of this cross section would be an extremely 
interesting and important probe of the Higgs self-couplings, especially given that all Higgs bosons other than the 
are too heavy to observe directly at the linear collider in the situations considered. 

In short, the lesson of this section is that if nature chooses the Higgs sector parameters to be far from the decoupling 
regime, it could happen that only the will be within the kinematic reach of the linear collider and that the vvJ^ 
cross section might be observable. If in the future the LHC finds a fairly heavy CP-even Higgs boson, then, within 
the 2HDM (or similar) context, the type of situation considered here will be required for consistency with current 
precision electroweak constraints and one should urgently search for the in single production modes. 



C. Special situations in the MSSM 

In order to obtain a large cross section for vv_fi^^ production in the MSSM when tan/3 > 1, the ti or 66 coupling 
of the A" must be enhanced very substantially. This is within the realm of possibility. In particular, at large tan/3 
it is possible to have important one loop modifications to the J^hh coupling due to radiative corrections involving a 
gluino and a bottom squark. We briefly review this aspect of the MSSM and then show the resulting enhancement 
for favorable parameter choices. 

Since supersymmetry is broken, the bottom quark will have, in addition to its usual tree-level coupling to the Higgs 
field <i>5i small one-loop-induced couphng to $2 that couples to up quarks at tree- level: 

- >CYukawa ^ h^^\hh + {^h^)'^\hh . (20) 

When the Higgs doublets acquire their vacuum expectation values, the bottom quark mass receives an extra contri- 
bution equal to (Ahb)v2- Although Ahi, is one-loop suppressed relative to hh, for sufficiently large values of tan/3 
{v2 ^ vi) the contribution to the bottom quark mass of both terms in Eq. (|20|l may be comparable in size. This 
induces a large modification in the tree-level relation, 

m, = ^(l + A,), (21) 

where Af, = (Ahb) tan/3//ib. The function At contains two main contributions: one from a bottom squark-gluino loop 
(which depends on the two bottom squark masses mg^ and mj^^ and the gluino mass nig) and another one from a top 
squark-higgsino loop (which depends on the two top squark masses m^^ and m^, and the higgsino mass parameter 
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FIG. 10: Left-hand graphs: a{uVA°) as a function of tan /3 for Ecm = 800 GeV, wfth (points slightly to the left) and without 
(points slightly to the right) Higgs-self-coupling triangle diagrams. Right-hand graphs: corresponding values of AT,(7/(47r). 
(We do not attempt to display point-to-point matching between the left-hand and right-hand plots.) The three plots are for 
m^o = 100, 350, 600 GeV, m^o = 800 GeV, rrijjo = 900 GeV and m^± G [900, 930] GeV. The plotted points have m„± and a 
chosen so that S, T lie within the 95% CL precision electroweak data ellipse. 
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FIG. 11: Dependence of a{uVA°) on tan/3 for m^o = 300 GeV and ^ = 800 GeV, using HDECAY and MSSM parameters as 
specified in the text. Results are shown keeping only the (dominant) f-channel (i.e., WW fusion) t and b fermion-loop graphs. 



/i) . The explicit form of Ab at one- loop in the limit of Ms ^ is given by |^ |^ : 



A,, 



2a, 
37r 



■mgfita,nf3 I{mi^,mi^_,mg) + — At^tan/J I{mf^,m^^, fi), 



(22) 



where = g'^/Air, Yf = h^/Ai:, and contributions proportional to the electroweak gauge couplings have been neglected. 
The function / is manifestly positive. Since the Higgs coupling proportional to A/if, is a manifestation of the broken 
supersymmetry in the low energy theory. A;, does not decouple in the limit of large supersymmetry breaking masses. 
Indeed, if all supersymmetry breaking mass parameters (and /i) are scaled by a common factor, the correction Af, 
remains constant. For our purposes, the important implication is the modified form of the bbA^^ coupling j/f, [compare 
Eq. ini)]: 



Vb 



rrib tan 13 
Vl + Ah 



(23) 



where A;, oc tan/3 [see Eq. (|22|) ]. For appropriate parameter choices with /i < (assuming the standard convention of 
TOg > 0), Ab ^ —1 will occur at some value of moderate to large tan/3. At and near this point, the triangle diagrams 
in our calculation involving the bbA^ coupHng [diagram (2) of Fig. ^and the 5-loop cases of diagrams (1) and (2) 
of Fig. E] will be greatly enhanced leading to a very large cross section. We illustrate this for the specific choices of 

ruAo = 300 GeV, ^ 800 GeV, n = -2 TcV, to- = 525 GeV, to- = M2 = to- = to- = to- = to- =1 TeV, 

and At = Ab = Ar = fi/ tan/3+ V6to~^ (corresponding to maximal-mixing in the stop sector). Since At 7^ in general, 
both the gluino-bottom squark and higgsino-top squark loops contribute to A;,. Note that we have chosen sufficiently 
large masses for the SUSY particles that the one-loop contributions to e+e^ vVA^ involving them will be very 
suppressed. Keeping only the WW fusion t and b loop diagrams [diagrams (l)-(2) of Fig. ^ for simplicity, we plot 
a{vT'A^) as a function of tan/3 for these parameter choices in Fig. and compare to the result that would be 
obtained without including A;,. A close-up of the tan/3 region in which A^ ^ — 1 is shown in Fig. E| The plots show 
that for tan/3 within a few per mil of the point where A5 ~ — 1, the cross section can approach the femtobarn level. 
However, this corresponds to an extremely nonperturbative bbA^ coupling ?;^/(47r) ~ 10^. Requiring perturbativity, 
?/^/(47r) < 0(1), yields a cross section of order lO^"' fb, an enhancement of 1-2 orders of magnitude compared to the 
cross section without the Ab effects. 
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FIG. 12: (Left) Dependence oi aiwA") on tan /3 as in Fig.[TTl for 24 < tan/3 < 26. (Right) The corresponding value of y^/(47r). 

While it would be rather serendipitous for the MSSM parameters to be within the rather narrow range of tan/3 for 
which Ah is sufficiently near —1 to yield a significantly enhanced e~''e~ vvA^^ cross section, one cannot simply rule 
the possibility out. Of course, if the one-loop enhancement is very large, higher loop corrections would need to be 
computed to more precisely evaluate the magnitude of the cross section enhancement. 



V. CONCLUSIONS 



Detection of the CP-odd of a CP-conserving two-Higgs-doublet model via tree-level production mechanisms 
might not be possible due to: (a) the absence of VVA^ {V = W, Z) tree-level couplings; (b) kinematic limitations 
such as ^/Sl,+e- < 2m^o,m^o -|- mHo; and/or (c) the small size of the "Yukawa radiation" processes yielding tiA^ 
and final states [as typical for intermediate tan/3 values in the "wedge" region of (mAo,tan/3) parameter space]. 
These difficulties become magnified in models with more than one CP-odd Higgs boson (the minimal such Higgs sector 
is that containing two-doublets plus one-singlet). Thus, it is generically important to compute single- A*^ production 
rates deriving from one-loop diagrams (i.e. not associated with A'^ radiation from a b or t quark) and ascertain the 
circumstances under which such processes might allow detection of the j4°. In this paper, we have computed the full 
(one-loop) cross section for e+e" vvAP in the general CP-conserving 2HDM. Our results are presented in such 
a way that they can be easily extended to more complicated models containing a CP-odd Higgs boson. Complete 
formulae for the matrix elements are given in the Appendices. 

We have presented numerical results for three cases in the context of the CP-conserving type-II 2HDM. The first case 
considered is that where the Higgs sector parameters are chosen using the tree-level MSSM Higgs sector constraints. 
For this choice, the 2HDM rapidly enters the "decoupling" regime once m^o > mz- For tan/3 > 1, the e+e~ — ^ vv^ 
cross section is typically rather small, especially for m^o > 2m(. However, for tan/3 < 1 we find that e+e~ vvAP 
production could provide a viable A" signal, thus covering this important part of the (TO^o,tan/3) "wedge" parameter 
space region in the 2HDM where discovery using tree-level processes is not possible. In addition, if detected the 
vvAP rate would provide a very sensitive constraint on tan /?. 

For the decoupling 2HDM parameter choices considered above, the diagrams contributing to vvJ^ production that 
involve Higgs self-couplings (A^A^ft," and A^A^if^) are typically smaller (often much smaller) in size than t- and 
6-loop diagrams involving tiAF' and couplings. Thus, we have explored alternative 2HDM parameter choices for 
which one is far from the decoupling limit and the Higgs self-couplings are as large as they can be without violating 
perturbativity or precision electroweak constraints. We have found that substantial enhancement of the e~ vVJ^ 
cross section is possible when tan/3 > 10, possibly sufficient to make the process marginally observable. While AP 
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discovery via the i/PA^ final state would probably still be difficult, if the has been detected by other means and the 
approximate value of tan/3 is already known, the above results imply that e~^e~ — > i/uA'^ might provide an especially 
sensitive probe of the Higgs self-couplings of the 2HDM. 

Finally, we considered the effect on the cross section of an enhancement of the bbA'^ coupling caused by SUSY 
radiative corrections in the MSSM, parameterized by A;,. While the At effects can enhance the cross section by 1-2 
orders of magnitude while maintaining a pcrturbative bbA'^ coupling, this enhancement typically occurs at moderate 
to large values of tan/3 where the cross section is already quite tiny, so that even with a A;, enhancement the cross 
section is not larger than about 10~^ fb. 

For most parameter choices, the e^e^ vvAl' cross section is smaller than other single A° production channels, 
most notably e~ — > 7^*^. However, the two production modes are complementary in important respects. First, the 
vvAl' cross section could provide confirmation of a signal seen in the 7^" final state. Second, if tan /3 is not large, as 
might be known either because the two rates are fairly large or from other Higgs sector measurements, then A;, or 
Higgs self-coupling enhancements cannot be substantial and the determinations of tan /3 provided by the two processes 
can be fruitfully combined. Third, while the 7^° rate could also be enhanced by Af, effects it cannot be enhanced by 
large Higgs self-couplings (the needed A^H^H^ and ZH^G^ vertices for the H~^-\oop and iJ^-G"*" -loop, respectively, 
being absent in the 2HDM). Thus, an unexpectedly large cross section in the e+e^ vvA^ channel would signal 
large Higgs self-couplings if a similar enhancement is not found in the ^A^ final state. We note that this cross-check 
would be important even if no evidence for SUSY particles is found since a large Af, can arise for arbitrarily large 
SUSY particle masses. 

An important extension of this work will be to include the contributions from one-loop diagrams that involve 
supersymmetric particles. For a light SUSY spectrum, very substantial enhancements could occur. This could be 
especially important in the following situation. Imagine that the LHC (or Tevatron) discovers fairly light SUSY 
particles and a SM-like /i° but is unable to detect the heavier Higgs bosons and H^. This situation arises 

at the LHC for moderate tan/3 values within the wedge beginning at m^o ~ 200 GeV and becoming rapidly wider 
in tan/3 as niAa increases. If a linear e^e~ collider has too low a center-of-mass energy for A^H^ and H~^H~ pair 
production, we must search for the A^ (and the and -ff*) in the single production modes. For a light enough 
SUSY spectrum these modes could be sufficiently enhanced to make 7^°, i^vA'^ and similar processes observable, as 
found to be the case for e+e~ W^H^ production [s^ . 
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APPENDIX A: NOTATION AND CONVENTIONS 

We follow the notation used on the LoopTools web page (as of the date of this paper) for the one-loop integrals. 
To avoid any possible confusion, our explicit conventions are given below. The two-point integrals are: 

lb - /|P - ^f)«,''+'^)' - .nit ■ <'^" 

where D is the number of dimensions. 
The three-point integrals are: 



16^2 ^ ' ' ' J {2t:)D {q^-ml){{q + kiY-ml){{q + k2Y-miy 

where the denominator structure follows from the Feynman diagram of Fig. 1131 The tensor integrals arc decomposed 
in terms of scalar components as 

fjt,. ^ gf^''Coo + Kk'(Cii + k^kl^C22 + {k'^k!^ + k^k'()Ci2. (A3) 
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FIG. 13: Illustration of the LoopTools conventions employed for the triangle diagram case. 



The arguments of the scalar three-point integrals are specified in our convention as (/c^, (^2 — fci)^, fc|, m^, m^, tri'^). 
The four-point integrals are: 

{Uq^V^.V^ I - / 77—7777-2 _2^//„ , ^2 _2^//„ , 7. \2 ^2^/^„ , ^2 ZriT' V^^) 



16^2 I ' J 7 (27r)^ (g2_„j2)((^^fc^)2_^2)((^^;.2)2_^2)((g^;,3)2_^2) 

where the tensor integrals are decomposed in terms of scalar components as 

jj^u ^ + fcf fc^^ll + ki^kl^D22 + k^k^D33 

+ (fcf fc^ + ki^k'i)Di2 + (k'^k^ + k^k-^Dis + (k^k^ + k^k^)D23. (A5) 

The arguments of the scalar four-point integrals are {kf, (^2 — fci)^, (^3 — ^2)^, fef, ^21 (^1 ~ fe)^7 "^ii '^12, m^, ml). 

For the three-point functions Ci and Cij, it is useful to define the sums and differences of one-loop integrals as 
follows: 

(^Lji^i^ (^2 - kif ,kl,rn\,ml,ml) = i [C,,y(fc^, (^2 - kif ,kl,m\,ml,ml) + Ci,y((/s2 - fci)^, fc^ , m^, mj, TO3)] 

(^2 - A:i)2,fc2, 7711,7712,7713) = i [Cj,y(fci, (fc2 - A:i)2,fc2, 7711,7772,7773) - Cj,y((/s2 - kif,kf,kl,ml,ml,ml)] 

(A6) 



APPENDIX B: 2HDM CONTRIBUTIONS TO T-CHANNEL DIAGRAMS 



We list here our results for the i-channel, i.e. VFW^-fusion, diagrams. In the expressions below, ki and fc2 denote 
the momenta of the and W~ , respectively, with directions such that fci = p,^ ~ Pe- ^^iid k2 = Pij — Pe+] see Fig. 
^1. (These fci_2 should not be confused with those defining the LoopTools conventions in Appendix 1X1 ) 

Fig. nil (tt6'loop): 



_ g3m[cot^ 



(7^771? cot /?iVr n 

= ^Stt ^^2"^, B3 



with the arguments for the integral functions as C{kl, k2,rn^o,rri^ ,rri^,rri^). 
Fig. 1112 {bbt loop): 



_ ^ g^ml tan pN, s 

- « — -z ^ — C2 , (B4) 

1077lVl/7r^ 
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g^mg tan ^TVc ^^ 
IGtovi/tt^ ^ ' 



^2)^*2 



(B5) 
(B6) 



with the arguments for the integral functions as C(fcJ, to^o, , 
Fig. [23 (S'5S' loop): 



Gsss 
Hsss 



27r2 



r5ASiS3ffH'S2Sifl'H'S3S2[C22 + C'12 + Q 



(B7) 
(B8) 



with the arguments for the integral functions as C(fcJ, m^o, "^2^ ^i)- 

The combinations of scalar particles to be summed over and the respective couplings are given for the general 
type-II 2HDM in Table Note that the "flipped" diagrams with 6*1 ^ 6*3 have already been taken into account in 



Si 


S2 


S3 


9A0SiSa 


9WS2S1 


gws^s^ 


h° 


H+ 


A° 


9 AO AO 


— cos(P — a) /2sw 


—i/2sw 


H° 


H+ 


A° 


9a<>aoh<> 


sin(/3 — a)/2sw 


—i/2sw 


h° 


G+ 


G" 


(m\o - mlo) cos{P - a)/2mwsw 


— sin(/3 — a)/2sw 


—i/2sw 


H° 


G+ 


G° 


("1^0 — 'rn^o) sin(/3 — a)/2mwsw 


— cos(/3 — a)/2sw 


—i/2sw 


H- 


ftO 


G- 


i{m'^j^± — ■m\a)/2mwsw 


cos(/3 — a)/2sw 


sin(/? — a)/2sw 


H- 


//" 


G- 


i{m'^jj± --m\o)/2mwsw 


— sin(/3 — a) /2sw 


cos(/3 — a)/2aw 



TABLE I: The A°SS and VSS couplings needed for Fig. 03 are tabulated. Here, sw and cw are the sine and cosine of the 
Weinberg angle. 



the form factors given above. 

In the MSSM, the and coupling coefficients are: 

dAOAOho = {—mz/2swcw)cos2l3B\n{l3 + a), 
9aoaoh° = {mz /2swcw) cos 2 fi cos{l3 + a). 



(B9) 
(BIO) 



In the general 2IIDM, these two coefficients are best expressed |2J| in terms of certain combinations of the A; of 
Eq. 1^. For Ae = A7 = 0, as assumed in this paper, we have 



9A°A°h° 



[ArS/3- 



9a°a°h° = --[Atc/3- 



A(7C/3-a] , 
Ac/S^-q] , 



where - 



2cw^z 
swg^ 



and 



At = is^^(Ai + A2) + A345(4 + 4) 

A;^^ = |s2/3(s|Ai — C^A2 + C2/3A345) . 



2A., 



(Bll) 
(B12) 

(B13) 
(B14) 



In Eqs. l)B13p and (|B14() the quartic Higgs couplings Ai^2,3,4 arc given in terms of the parameter set m^o, m/fo, mho, 
mH±, a, tan/3 and A5 by Eqs. of Scc.lTll 

Fig.[Tl4+Fig.[Tl5 {SSV loop): 



Gssv 



-^9ASiV39WS2SigwV3S2^00^ 



Hssv — -r^ 9 AS1V39W S2Si9wv3S2[C22 + - C2 ] 



An 



(B15) 
(B16) 



with the arguments for the integral functions as C{ki,k2,m\o,m\,m^,'m^). The combinations of scalar and vector 
particles to be summed over and the respective couplings are given in Table ^1 Again wc have already taken into 
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Si S2 Vs 




9WS2S1 


gwv3S2 


H- h° W- 
R- H° W- 
ftO G+ Z 
H° G+ Z 


i/2sw 
i/2sw 
[—i/2swcw) cos(/? — a) 
{i/2swcw) sm(/3 — a) 


cos(/3 — a)/2sw 

— sin(/3 — a)/2sw 

— sin(/3 — a) /2sw 

— cos(/? — a)/2sw 


(mw/sw) sin(/9 — a) 
(mw/sw) cos(/3 — a) 

—mzsw 

—mzsw 



TABLE II: The VSS couplings needed for Figs.[T]4 and[T]5 are tabulated. 



account the "flipped" diagrams with S*! <-> V3. 
Fig.[Tl6+Fig.[Il7 {SVV loop): 



Gsvv = — SASiVaSwyaSi 1 [Cqo + • ^2^1 +m\aC2 -m\C^ 



Hsvv 



le cw 
Stt^ Sw 



9AS^V39wv2St [C22 + - 4Cf - C2] 



(B17) 
(B18) 



with the arguments for the integral functions as C(fc^, fc|, m^o , ti^, TO2, m§). The combinations of scalar and vector 
particles to be summed over and the respective couplings are given in Table ITTTl Again we have already taken into 



5i V2 V3 


9 AOs-, 


9WV2S1 


h" W+ Z 
H° W+ Z 


{—i/2swcw) cos(/3 — a) 
{i/2swcw) sm{P — a) 


(mw/sw) sin(/3 — a) 
(mw/sw) cos(/3 — a) 



TABLE III: The A°SV and WVS couplings needed for Figs.[T]6 and[Tl7 are tabulated. 

account the "flipped" diagrams with S*! <-> V3. 

It is easy to check that in the on-shell limit where kl — k"^ — rn^, all the C/^^^ 's go to zero. Contributions to G 
and H vanish and the only contribution to W'^W~ effective coupling comes from F, as pointed out in [i^. 

Fig.njS (box diagrams): 



G 



Hi — H4 



9 

Y^9z9wwsgzsA{k2 - m'^-) [Cq + £'o(m| + 2s + 2t2 + 2ui) 

+Di{2m\o +s + t2 + ui) + D2{2m\o -2s- 2ti - 2u2 + 12) + A3(2s + 2t2 + 2ui)] (B19) 



9" 



-.9'z9wwsgzsA{k\ - m'^){4D2), 



(B20) 



with the arguments for the integral functions as I?(to^o, ^i, 0, 0, ^2, s + ti + it2, m^, m^, 0), Co(ii,0,s + ti + 
M2,m.|,"i^,0), gz = -T3^/cw, and 

S = 2pe- • Pe+ , S^ 2py • pu, ti^kl^ -2pe- • p^, t2 = kj = -2pe+ ' Pu, = -2pe- ' PF, "2 = -2pe+ • Pty ■ 

(B21) 



The scalar could be h'^ or H'^, with the couplings given in Table ITVI 



s 


9WWS 


9ZSA0 




mw sin(/3 — a) 
mw cos(/3 — a) 


cos(/3 — a)/{2cw) 
-sm{l3-a)/(2cw) 



TABLE IV: The WWS and ZSA° couplings needed for Figs. [T]8-[T]ll are given. 

Fig.n]9: Similar to Fig^S, under the exchange of 

H4 H^, ui ^ U2, ti ^ t2, k\ ^ fc|, and an overall " — " sign. 
Fig.QJlO: Similar to Fig^S, under the exchange of 

9z g^z === -~{T-ieL - QeSw)/(^w, Hi H2, H4 ^ H3, s ^ s, ui ^ U2, and an ovcraU 



(B22) 
sign. (B23) 



21 



Fig.njll: Similar to Fig. ^8, under the exchange of 



9z 



eL 

9z 



-{TaeL - QeSw)/cw, Hi ^ H2, s ^ s, ti^t2, kl^kl 



(B24) 



APPENDIX C: 2HDM CONTRIBUTIONS TO S-CHANNEL DIAGRAMS 



For the s-channel diagrams, we introduce the foUowing operators: 



O3 ^V{p^+) P^Pru{p^-)u{p^) ^e-PLv{pj7) 

C5 ^v{Pe+) i)yPRu(jp^-)u{p^) ^e+PLv{pj7) 

07^v{Pe+) ^PRu{p^-)u{pt,) ^e-PLv{pj7) 

O9 ^~'{Pe+) i>i7PRu{p^-)u{pi,) P^+Plv{pj7) 

We now hst our results for s-channel diagrams. 
Fig.[21l4|ll2 (s-channel top quark loop): 



O2 = v{Pe+)l^PLu{pe-)u{pu)lf,PLv{pT7) 
O4 ^v{p^+) i)yPLu{pe-)u{Pv) A'-Plv{p17) 

Ce ^v{p^+) i)vPLu{pe-)y-{Pv) A+Plv{pj7) 

Os =v{Pe+) i>i7PLu{pe-)u{py) A- Plv{p17) 
OlO =v{pe+) i>vPLu{Pe-)u{pu) i>e+PLv{pu)- 



(CI) 



M = 



em^ cot (3 



w 

(ti -t2+ui- u2){-gf Oi + 5^02) - gv^'iOs + 0r~05- O9) + .9^^(04 + Os - Oe ~ 



{s — ■my){s — m^^ + imzT z) mws 
1 
2 



no) 



,(C2) 



with the arguments for the integral functions as C{m?^Q, s, s, , , ml). This includes the top quark going clockwise 
and counterclockwise. The gauge boson connecting the initial e+e^ to the top quark loop is 1/ = 7 or The couplings 
are defined as: 



qL 

9z 



qR 

9z 



Qq^W 



9t 



To get the s-channel bottom quark loop, make the following substitutions: 

mt^rrib, cot/3 ^ tan/3, (# ,5^^) {gv.9f)- 
Fig. [33 (box diagram): For /i" in the box. 



-Qq- 



(C3) 
(C4) 



M 



o? sin 2(/3 — a)mw 



e9z 



aJ2 ^3 



X [iiaffOi + {gffO^) (-Co + {s + h+ui- mlo)Do 



(s — m| + imzTz) 

+ {h + ui)Di + {m%, -s-h- ui)D3) + ADiiigffiO^ + Or) + {gf)''{04 + Og))] , 



(C5) 



with the arguments for the integral functions as C(0, 0, s, m^, 0, m?z) ^^'^ -D(s, 0, 0, ?ti^07 (s+^i+iti), s, m\o-: 'ti^, 0, m\). 
The crossed box is obtained by applying the substitutions: 



ti^t2, 1*1^^2, (03,04,07,08)^(0.5,06,09,010), and an Overall "-" sign. 
The box diagram containing H'^ is obtained by applying the substitution: 

m/jo — > muo and an overall "— " sign. 

Fig. 124: For /i° in the box, 

sv[i2{(i - a)mw e{gz^)^ 



(C6) 
(C7) 



M 



(s-m|) 



4-(t2 +Ml)£'l + ("1^0 



[(.gf Oi + gf02){-C^ + (s + + ^.ti - ™'o)i?o 
t2 - ^i)i?3) + 4i?i(.gf (O7 + O9) + #(08 + Oio)] 



(C8) 



with the arguments for the integral functions as C(0, 0, s, m^, 0, ?n^) and I?(s, 0, 0, m^o, u', s, ti^o, to^, 0, m^), where 
m' = s -I- 12 + ui . 

The crossed box is obtained by applying the substitutions: 

t2^ii, ui^U2, (O7,08,09,0io) ^ (03,04,05,06), and an overall "-" sign. (C9) 

The box diagram containing is obtained by applying the substitution: 

m/jo — ^ m/fo and an overall "— " sign. (CIO) 
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APPENDIX D: SQUARE OF THE MATRIX ELEMENT 

The cross section for e+e^ vvAl' is evaluated by integrating the spin-averaged matrix element square [Eq. (|DT|| ] 
over the three body phase space of the final states: 

lE . i-^P - . 1-^*1' + !e . i>'^(^z4)P + |e . 

— ^ spins 4:^ — spins 4:^ — ^ spins 4^ — ^ spins 



-E 



MtMsie-^elY +h.c., (Dl) 



where the 3 in the second and third terms represents the sum over three neutrino flavors for the s-channel contribution. 
The various pieces in Eq. IjDip arc given explicitly below. 

The spin-summed amplitude squared for i-channcl diagrams is 



E - 1-^*1' = I 

' ^ spins 4 



.9^ 



2 

Kt (D2) 



where 



+ i\G\^UiU2 

+ \Hi\'^s'^UiU2 + \H2\'^s^uiU2 + iH^l'^ssuj + |i74pSSU2 

+ 2'iR{HiH;)ssuiU2 

+ di{HiH^)sui{s's + U1U2 — tit2) — 43(iJiiJ3 )e(e~, e"*", 77, v)sui 

+ ^{HiHl)su2{ss + uiU2 - tit2) + ^'^{HiHl)e{e-,e+,V,iy)su2 

+ ^{H2H^)sui{ss + U1U2 - ii<2) - 4Q?(iJ2-ff3)e(e", e+, i7, v)sui 

+ 'iR{H2Hl)su2iss + U1U2 - tit2) + 43(i72i?4*)e(e-, e+, j7, v)su2 

+ n{H3H^){s^s^ - 2sstit2+tltl - 2tit2UiU2 + ulul) + ^'^{H3Hl)e{e- ,e+ ,V,v){ss + U1U2 - ^1^2) 

+ 3fi(i^G*)8e(e",e+,F,z/)(wi +U2) 

-|- 23(_FG'*)(— SSUi + SSU2 + tit2Ui — t\t2U2 + u\u2 — Uiu\) 

-\- 4M{FHi)e{e'~ , e^,F, j/)(sui -f- SU2) + s^sui + s^su2 -t- su^W2 — suim| -|- stit2Ui — stit2U2) 

+ 4M{FH2)e{e~ , ,V, iy){'sui + SU2) + '<s{FH2){—ss^ui + ss^U2 + ^u\u2 — 'suiu\ + siii2Mi — siit2U2) 

+ 43fJ(Fi?*)e(e-,e+,I7,j/)(s?+u? -tit2) 

+ '^{F H'^)[u\u2 — tit2u\ — SSu\ — tit2UiU2 + t\t\ + S^J^ — 2sstit2 — SSU1U2) 

+ 45R(FiJ*)e(e-,e+,I7,i/)(s?-f 

— ^{F H1)[uiu\ — tit2u\ — SSU2 — tit2U\U2 + iii2 + ~ 2sstit2 — s'suiU2) 

+ A^{GHl)suiU2 + m{GH;)suiU2 

+ 2^{GH;)ui{uiU2 - ht2 + ss) - 85(G'i?3)e(e~, e+, j7, v)ui 

+ 23fi(GiJ4)u2(wiW2-ii<2 + s?) + 83(Gi?4)e(e~,e+,i7,iy)u2, (D3) 



and 



with e0i23 _ Here, and 5 denote the real and imaginary parts, respectively, of the indicated products. 
The spin-summed amplitudes-squared for s-channel diagrams are 

\Ms{e^e+)\^ = A\Mi\^ht2 + \Ms.\HluiU2 + \M5\^tit2ul + \M7\^ht2ul + \Mg\hluiU2 



spins 



+ 2^{MiMl)ti{ss~tit2-uiU2)-^'^{MiM%,)tie{e-,e+,V,v) 

- A^{MiMl)ht2U2 - m{MiM^)ht2Ui 

+ 23?(Mi7W;)t2(ss - tit2 - uiU2) + %'^{MiMl)t2e{e- ,e+ ,V, v) 

- '^{MzMl)tiU2{ss - tit2 - U1U2) - 43(7W3A^5)iiU2e(e",e+,F,i/) 



'Si{M3MY)tiUi{ss - tit2 - U1U2) - 43(A^3>ly)tiWie(e ,e+,i7, j/) 



-{ss ~ tit2 - uiu^f - 8(e(e , e+, v, v)f 



+ 43(7^3X9) (ss - <it2 - MiM2)e(e", e+,T7, j/) 



-(s.s - tit2 - MiU2)^ + 8(e(e-, e+, 17, v)f 



'^{M^Ml)t2U2{ss - tit2 - U1U2) - 43(A^5A^;)t2W2e(e ,e+,i7,i/) 
3fi(7W7Ai;)i2Wi(ss - tit2 - U1U2) - 43(7W7A^9)t2Uie(e", e+,77, ly) 



4|A^2pwiW2 



|A^4ptlltlU2 



|Al6|'ilt2U2 



- 45R(7W2Al4)^i'^iW2 

+ 2n{M2Ml)u2{ss ~ tit2 - uiU2) + 83(7W2X;)u2e(e-, e+,I7, i/) 
+ 2'^{M2Ml)ui[ss - tit2 - U1U2) - %'^{M2Ml)uie{e~ ,e'^ ,v,u) 



m{M2M*w)t2UiU2 
^{M4,Ml)tiU2{sS - tit2 

'?R{M4M8)tiui{ss - tit2 



-{ss - tit2 



77 (ss 



U1U2) — 4Q(A^4A^g)tiM2e(e ,e^,v,v) 
U1U2) +45(A^4A^8)tiUie(e~,e+,F, v) 



U1U2) + 8(e(e ,e+,t/, i/)) 



U1U2Y ~ 8(e(e" 



A'^{MqM*8){ss — <it2 — MiU2)e(e , e+, i/, i^) 
5R(7W6Xto)^2W2(ss - tit2 - U1U2) + 43(7W6Xto)i2M2e(e" 
5R(7W8Xto)i2Ui (ss - iits - U1U2) - 43(A^8A^to)*2Mie(e", e+,F, i^), 



where 8(e(e , e+, i^, i^))^ ~ ^\{^^ + Wiit2 — tit2)^ + 2s.suiW2. 
The interference terms between s- and t-channel diagrams are 



V MtMs{ere+)* +\\.c. 

^ — ^ spins 



I even spins 
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2 (ti-J7l2^)(t2 



where 



~m{GM*2)uiU2 - m{FM*2){ui+U2)e{e- ,e+,V, u) 
-23(FX;)(ui - it2)(-ss + tit2 + U1U2) - A^{HiM*2)suiU2 
-4,^{H2M*2)suiU2 - 2'^{H3M*2)ui[uiU2 - ^1*2 + ss) 
-83(i73Ai;)uie(e-, e+,77, i/) - 2^{H iM*2)u2{uiU2 - iit2 + ss) 
+8'^{HAM*2)u2e{e- ,V, u) 
A^{FMl)h{ui + W2)e(e-, e+,T7, u) 

-Q{FM4)ti{sS - tit2 - UiU2){ui - U2) 

+45R(GA^*)tiMiW2 + 2^{HiMl)stiuiU2 + 2^{H2Ml)stiuiU2 
+^{H3Ml)tiUi{ss + U1U2 - ht2) + 4':s{H3Ml)tiUie{e-,e+,V, ly) 
+^{H4Ml)tiU2iss + U1U2 - tit2) - 43(i74A<4)iiU2e(e", e+,77, jy) 
45R(i^7Wg)(tit2 + U1U2 - ss)e(e", e+,77, jy) 
+'^{FM*e)[{ht2 + U1U2 - ssf - 2tit2U2{ui + U2)] 
+2^{GMl)u2{tit2 + U1U2 - ss) - d,'^{GMl)u2e{e- ,e+,V, v) 
+^{HiMl)su2{tit2 + U1U2 - ss) ~ 45(iJiAl6)su2e(e", e+,I7, i^) 
+^{H2Ml)su2itit2 + U1U2 - ss) - 43(iJ2A^6)'5"2e(e", e+,77, u) 

+^{H3Ml)[-iht2 - SS)2 + UiU2{tit2 + ss)] 
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+^iHiMl)uliuiU2 - ht2 - ss) - 4^{HiM*e)ule{e-,e+,V, u) 
Kg = m{FM*s){ht2 + uiU2- ss)e{e-,e+,V,u) 

+^{FMl)[-itit2 + UiU2 ~ Ssf + 2tit2Ui{ui + U2)] 

+2^{GMl)ui{tit2 + U1U2 - ss) + 83(G'7W;)wie(e-, e+, F, v) 
+3i{HiMl)sui{tit2 + U1U2 - ss) + 43(iJi7Wg)sMie(e", e+, 77, v) 
+^{H2Ml)sui{tit2 + U1U2 - ss) + 45(iJ2A4g)suie(e", e+,I7, v) 
+3i{H3Ml)ul{uiU2 - tit2 - ss) + 43(iJ3A^J)w^e(e~, e+, j7, v) 
+^{HiMl)[-{tit2 - ss)^ + uiU2{tit2 + ss)] 
+4^{H4M;){ss - tit2)e{e-,e+,V, v) 
K[a = A^{FM\Q)t2{ui + U2)e{e-,e+,V,v) 

~'^{FM*iQ)t2{sS - tit2 - UiU2){ui - U2) 

+m{GMlQ)t2UiU2 + TR{HiMlQ)st2UiU2 + 2^{H2MlQ)st2UiU2 

+^{H-iMlo)t2Ui{ss + U1U2 - ht2) + 43(i73A^Jo)t2Mie(e-, e+, 17, v) 
+n{HiMlo)t2U2{ss + uiU2-ht2)-4'^iHiMlo)t2U2eie',e+,V,iy). (D8) 



J. F. Gunion, H. E. Haber, G. L. Kane and S. Dawson, The Higgs Hunter's Guide, (Addison- Wesley, Redwood City, CA, 
1990) [Erratum arXiv:hep-ph/9302272 . 

LEP Higgs Working Group, LHWG/2001-03 (July 2001), arXiv:hep-ex/0107029 
LEP Higgs Working Group, LHWG/20Q1-Q4 (July 2001), arXiv:hep-ex/0107030 

For a recent summary, see e.g., U. Schwickerath, (Final) Higgs Results From iiSP arXiv:hep-ph/0205126 
For a recent summary, see e.g., L. Moneta, Higgs Searches at the Tevatron, to appear in the proceedings of 36th Rencontres 
de Moriond on QCD and Hadronic Interactions, Les Arcs, France, 17-24 Mar 2001 arXiv:hep-ex/0106050 . 
Y. Yamada, K. Hagiwara, and S. Matsumoto, Prog. Theor. Phys. Suppl. 123, 195 (1996) arXiv:hep-ph/9 512227| ; J. Erler 
and D.M. Pierce, Nucl. Phys. B526, 53 (1998) arXiv:hep-ph/9801238 . 

C. A. Marin and B. Hoeneisen, hep-ph/0210167 

B. Grzadkowski, J. F. Gunion and J. Kalinowski, Phys. Rev. D 60, 075011 (1999) |arXiv:hep-ph/9902308| . 
OPAL CoUaboration, G. Abbiendi et al, OPAL Physics jfete FN475 (2001), 
http : //opal .web. cern. ch/0pal/pubs/physnote/html/pn475 .html. 

M. Carena et. al, Report of the Tevatron Higgs working group arXiv:hep-ph/0010338 . 

ATLAS Collaboration, Detector and Physics Performance technical Design Report Vol. H (1999), CERN/LHCC/99-15, 
p. 675 - 811, available from http://atlasinfo.cern.ch/Atlas/GROUPS/PHYSICS/TDR/access.html. 

CMS Collaboration, Technical Design Report, CMS TDR 1-5 (1997,1998); S. Abdullin et al. [CMS Collaboration], Discovery 
potential for supersymmetry in CMS, J. Phys. G 28, 469 (2002) ,arXiv:hep-ph /9806366| . 

K. Lassila-Perini, ETH DisseHation thesis No. 12961 (1998). 

D. Zeppenfeld, R. Kinnunen, A. Nikitenko and E. Richter-Was, Phys. Rev. D 62, 013009 (2000) "arXiv:hep-ph/0002036" . 
M. Battaglia and K. Desch, in Physics and experiments with future linear e^e~ colliders, Proc. of the 5th Int. Linear 
Collider Workshop, Batavia, Illinois, USA, 2000, edited by A. Para and H. E. Fisk (American Institute of Physics, New 
York, 2001), pp. 163-182 ,arXiv :hep-ph/0101165 . 

J. A. Aguilar-Saavedra et oTpCFA/DESY LC Physics Working Group] , TEST A Technical Design Report, Part 3: Physics 
at an e"''e~ linear collider 'arXiv:hep-ph/0106315 . 

T. Abe et al. [American Linear Collider Working Group C ollaboration], L inear collider physics resource book for Snowmass 
2001, Part 2: Higgs and supersymmetry studies, arXiv:hep-ex/0106056 . 
J. R. Espinosa and J. F. Gunion, Phys. Rev. Lett. 82, 1084 (1999) arXiv:hep-ph/9807275 . 

G. W. Bennett et al. [Muon g-2 Collaboration], Phys. Rev. Lett. 89, 101804 (2002) JErratum-ibid. 89, 129903 (2002)] 
arXiv:hep-ex/0208001 . 

[20] T. Ibrahim and P. Nath, Phys. Rev. D 58, 111301 (1998) [Erratum-ibid. D 60, 099902 (1999)] |arXiv:hep-ph/9807501| ; 
M. Brhlik, G. J. Good and G. L. Kane, Phys. Rev. D 59, 115004 (1999) arXiv:hep-ph/9810457 ; M. Brhlik, L. L. Everett, 
G. L. Kane and J. L ykken, Phys. Rev. Lett. 83, 2124 (1999) arXiv:hep-ph/9905215 ; T. Ibrahim and P. Nath, Phys. Rev. 
D 61, 093004 (2000) |arXiv:hep-ph/99'T0553| ; T. Ibrahim, U. Chattopadhyay and P. Nath, Phys. Rev. D 64, 016010 (2001) 

trXiv:hep-ph/0102324r 
. E. Habcr and Y. Nir, Phys. Lett. B 306, 327 (1993) |arXiv: hep-ph/9302228 . 
[22] T. Farris, J. F. Gunion and H. E. Logan, in Proc. of the A PS/ DPF/ DPB Summer Study on the Future of Particle Physics 
(Snowmass 2001 ) ed. R. Davidson and C. Quigg, arXiv:hep-ph/0202087 



25 



arXiv:hep-ph/0110289l 



arXiv:hep-ph/9902417] 



[23] A. Djouadi, P. M. Zerwas and H. E. Haber, Multiple production of MSSM neutral Higgs bosons at high-energy e+ e- 
colUders, in Physics with e+ e- Linear Colliders (The European Working Groups 4 Feb - 1 Sep 1995: Session 1), Annecy, 
France, 4 Feb 1995, pp. 89-103, arXiv:hep-ph/9605437 

J. F. Gunion and H. E. Haber, Phys. Rev. D 67, 075019 (2003) •arXiv:hep-ph/0207010'. 

P. Chankowski, T. Farris, B. Grzadkowski, J. F. Gunion, J. Kalinowski and M. Krawczyk, Phys. Lett. B 496, 195 (2000) 
arXiv:hep-ph/0009271 . 

K. Cheung, C. H. Chou and O. C. Kong, Phys. Rev. D 64, 111301 (2001) arXiv:h ep-ph /0103183| . 

A. Djouadi, J. Kahnowski and P. M. Zerwas, Mod. Phys. Lett. A 7, 1765 (1992). 

B. Grzadkowski, J. F. Gunion and J. Kahnowski, Phys. Lett. B480, 287 (2000) hep-ph/0 00T093l. 
A. Gutierrez-Rodriguez, M. A. Hernandez-Ruiz and O. A. Sampayo, Rev. Mex. Fis. 48, 413 (2002) 
U. Cotti, A. Gutierrez-Rodriguez, A. Rosado and O. A. Sampayo, Phys. Rev. D 59, 095011 (1999) 

J. F. Gunion and H. E. Haber, Phys. Rev. D 48, 5109 (1993). 

D. M. Asner, J. B. Gronberg and J. F. Gunion, Phys. Rev. D 67, 035009 (2003) arXiv:hep-ph/01 10320 . 
M. M. MuhUeitner, M. Kramer, M. Spira and P. M. Zerwas, Phys. Lett. B 508, 311 (2001) arXiv:hep-ph/01010W . 
M. Carena, H. E. Haber, H. E. Logan and S. Mrenna, Phys. Rev. D 65, 055005 (2002) [Erratum-ibid. D 65, 099902 (2002)] 

arXiv:hep-ph/0106116 . 

M. Battagha and K. Desc h,|arXiv:hep-ph/0101165| 
S. Kiyoura and Y. Okada, arXiv:hep-ph/0101172" 

K. Abe et al. [AGFA Linear CoUider Working Group Collaboration], "Particle physics experiments at JLC," 

arXiv:hep-ph/0109166 

A. Djouadi, V. Driesen, W. Hollik and J. Rosiek, NucL Phys. B 491, 68 (1997) |arXiv:hep-ph/9609420| . 
A. G. Akcroyd, A. Arhrib and M. Capdequi Peyranere, Mod. Phys. Lett. A 14, 2093 (1999) [Erratum-ibid. A 17, 373 
(2002)] arXiv:hep-ph/9907542 . 

A. G. Akcroyd, A. Arhrib and M. Capdequi Peyranere, Phys. Rev. D 64, 075007 (2001) [Erratum-ibid. D 65, 099903 
(2002)]. 

A. Arhrib, Phys. Rev. D 67, 015003 (2003) arXiv:hep-ph/0207330^. 

T. Farris, J.F. Gunion, H.E. Logan, S. Su, in preparation. 

P. Chankowski, S. Pokorski and J. Rosiek, NucL Phys. B 423, 437 (1994) |arXiv:hep-ph/9303309| . 
J. F. Gunion, H. E. Haber and C. Kao, Phys. Rev. D 46, 2907 (1992). 

A. Denner, S. Dittmaier, M. Roth and M. M. Weber, Nucl. Phys. B 660, 289 (2003) arXiv:hep-ph/0302198". 
T. Hahn and M. Perez- Victoria, Comput. Phys. Commun. 118, 153 (1999) arXiv:hep-ph/9807565 ; T. Hahn, LoopTools 
User's Guide, http://www.feynarts.de/looptools/. We note that the notation in the article is not always the same as 
that which appeared on the web page at the time this paper was written. We have employed the web-page notation as 

explicitly given in the Appendix A text. 

S. Heinemeyer, W. Hollik and G. Weiglein, Eur. Phys. J. C 9, 343 (1999) 'arXiv:hep-ph/98 12472| . 
A. G. Akeroyd and A. Arhrib, Phys. Rev. D 64, 095018 (2001) arXiv:hep-ph/0107040 . 

R. Hempfling, Phys. Rev. D49, 6168 (1994); L. Hall, R. Rattazzi and U. Sarid, Phys. Rev. D50, 7048 (1994) 

hep-ph/9306309 . 

M. Carena, M. Olechowski, S. Pokorski and C.E.M. Wagner, Nucl. Phys. B426, 269 (1994) |hep-ph/9402253| . 
D. Pierce, J. Bagger, K. Matchev, and R. Zhang, Nucl. Phys. B491, 3 (1997) hep-ph/96062TiT 

H. E. Logan and S. Su, Phys. Rev. D 66, 035001 (2002) arXiv:hep-ph/0203270 ; Phys. Rev. D 67, 017703 (2003) 
|arXiv:hep-p h/0206135|. 



